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Abstract. The authors apply the generalized master equation to analyze time- 
dependent transport through a finite quantum wire with an embedded subsystem. 
The parabolic quantum wire and the leads with several subbands are described by 
a continuous model. We use an approach originally developed for a tight-binding 
description selecting the relevant states for transport around the bias-window defined 
around the values of the chemical potential in the left and right leads in order to capture 
the effects of the nontrivial geometry of the system in the transport. We observe a 
partial current reflection as a manifestation of a quasi-bound state in an embedded 
well and the formation of a resonance state between an off-set potential hill and the 
boundary of the system. 
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1. Introduction 

Studying the dynamics of a few-level quantum system in contact to its environment is an 
old problem in statistical mechanics and solid state physics. The typical example comes 
from quantum optics where an atom is exposed to electromagnetic radiation described 
by a quantum or semiclassical bosonic bath (reservoir) [Ij. The natural theoretical tool 
for investigating transitions and computing life times is the reduced density operator 
(RDO) method which provides dynamical information about the atomic system while 
averaging out the reservoir degrees of freedom [2j . The time-evolution of the occupation 
probabilities is given by the so-called master equation. When writing the equation of 
motion for the off-diagonal matrix elements of the density operator one ends up with a 
generalized form of the master equation whose mathematical structure was investigated 
by many authors, (see e.g. [3ll4j and the references therein). 

As the system-reservoir picture is rather general the same method can be used to 
investigate electronic transport at the nanoscale: one considers a mesoscopic structure 
which opens to particle reservoirs (leads) and looks for the electronic flow through the 
nanosystem in the presence of a bias and/or time-dependent flelds. In this setup the 
contacts between the leads and the sample play the role of the system-reservoir coupling 
Hamiltonian. Since the reduced density operator method focuses on the dynamics of the 
sample it can be used to describe the transient regime, pump-and-probe experiments 
[5j and counting statistics [6] . Theoretical calculations for few-level quantum dots have 
been performed by several authors [3 El E]. Bruder and SchoUer [H] solved a 
quantum master equation for the diagonal elements of the density matrix, while Konig et 
al. [12] developed a real-time diagrammatic technique for the reduced density operator. 
Usually one assumes that the memory eflFects can be neglected and looks for steady-state 
currents within the Born-Markov and rotating wave approximation (RWA). In particular 
in the RWA the diagonal and oflF-diagonal elements of the reduced density operator are 
decoupled [7]. Moreover, in the high bias limit the generalized master equation reacheas 
the simpler form obtained previously by Gurvitz Non-Markovian effects were also 
considered in more recent works [151 US [El- In particular Vaz et al. fT6] took 
advantage of the Laplace transform and wrote down the Redfleld tensor for a two-level 
system. Obviously, such a task becomes too difficult for a more complex systems. A 
scheme including interaction effects was recently presented by Welack et al. [I8j . Various 
ways to solve the generalized master equation were discussed and compared in the recent 
paper of Timm [19j . 

In the references mentioned above one deals with rather simple two-level systems 
and their spectral properties do not follow from a specific geometry. Also, for very 
small quantum dots the precise location of the contacts is not important so it is 
reasonable to consider a tunneling Hamiltonian that does not contain information about 
the localization of the coupled states. In our recent paper [20j we have solved a non- 
Markovian generalized master equation (GME) for many-level systems described within 
a lattice model, paying special attention to the geometry of the system. More precisely. 
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under a strong perpendicular magnetic field the energy spectrum of the system is a 
Hofstadter spectrum with edge and bulk states. Given the very different nature of 
the sample states one expects them to carry different currents and (in the transient 
regime) even with a clear dependence on the contact region. We have therefore explicitly 
constructed a tunneling Hamiltonian that takes care of the location of the contacts and 
also depends on energy. We have analyzed transient currents and their dependence on 
various parameters of the system as well as on the initial configuration. 

The aim of this paper is to implement the same method for continuous models and 
to focus on the internal dynamics of the system. As an application we consider both 
a pure finite quantum wire with parabolic confinement and a wire with an embedded 
subsystem (a Gaussian well or potential barrier). The finite wire is connected to semi- 
infinite leads of the same width. We attempt to identify effects due to the underlying 
subband structure and also the formation of bound states due to the presence of the 
embedded potentials. Another motivation of this work is to compare the results of the 
present method with the ones obtained previously via the time-dependent Lippmann- 
Schwinger formalism [211 [22] . Although one expects serious technical problems in the 
continuous model due to the large number of states and the quite complex form of the 
tunneling term we show here that one can actually say a lot about the time-dependent 
transport in extended systems by selecting a set of single particle states that are expected 
to be relevant for the transport process. 

The paper is organized as follows: Section II presents the model and the main 
equations, Section III contains the numerical simulations and their discussions while 
Section IV summarises the results. Some technical details are left to Appendix A. 

2. Model 

We consider a two-dimensional wire in the x^-plane. In the ^-direction the electrons in 
it are parabolically confined with the characteristic energy HQq^ but in the x-direction 
they are confined by hard walls at x = ±Lx/2. The corresponding single-electron 
Hamiltonian is 

hs = ho + V=^ + ^m^ny + y(r), (1) 

where V{r) is a potential representing a subsystem embedded in the wire. The 
eigenfunctions of hs are denoted by ^^^l^)' with a referring to the two quantum numbers 
of the system, n^, and Uy. The eigenvalues of hs are denoted by Ea- Semi-infinite leads of 
same parabolic confinement are attached to the finite wire or the system at x = ±Lx/2. 
The single-electron Hamiltonian of the left and the right leads is noted by /il or /ir, 
respectively. Their eigenfunctions are ^^^'^(r) and the eigenvalues are e^'^(g). Due to 
the subband structure of the leads the quantum number q stands both for the continuous 
wavenumber and the subband index n^'^. The semi-infinite leads have a hard wall at 
X = ±La^/2, but at time t = to the system is opened by coupling it to the leads in a 
time dependent fashion described below. Further information about the single electron 
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states of the quantum wire and the leads is in Appendix A. The parabohc confinement 
of the system and the leads brings with it a natural lengthscale with a^f^o = h/rrf . 

Using the eigenfunctions introduced above one can write the Hamiltonians of the 
disconnected subsystems in the spectral representation 

a q 

In order to describe the coupling between the two subsystems we shall add an oS- 
diagonal pertubation to hs + hL + 

l=L,R a q 

where the coefficients T^^ are meant to describe a coupling between pairs of eigenstates 
{^^^, ^^^} and will be defined below. The time-dependent part of the coupling is regulated 
by a switching function X^'^(^)- 

In light of the variable number of electrons in the open system as a function of time 
and their correlations caused by the coupling to the leads we resort to a many-electron 
description introducing creation (annihilation) operators for electrons in the leads c^^ 
(cqi) and in the system d'l {da)^ with / = L,R. The many-electron Hamiltonian of the 
system and the leads is then 

H{t) = Y,Ejj,+ J2 e\q)c\,c,i + H^{t), (4) 

a q,l=L,K 

where the standard tunneling Hamiltonian H^{t) = HJ^{t)+H^{t) describes the coupling 
of the system to the left and right leads 

H^it) = X\t) J2 {^Uida + {TiaTdlc^i] . (5) 

In our previous work Ref. ([20j) where the system is described on a lattice we have 
explicitly constructed the coefficients T^^ for a neareast neighbor coupling appropriate 
for the lattice model. More precisely, by denoting by the site of the one-dimensional 
lead / which is to be coupled to its neighbor site ii of the sample T^^ was introduced as 
follows 

(X {^[{oW'^li^i)- (6) 

Let us now construct an appropriate generalization of Eq. ^ in the context of the 
continuous model under study. Basically, since we need the states of the subsystems to 
be mutually orthogonal in order to derive the GME we shall integrate a non-local overlap 
of a pair of eigenstates (^^^, ^^^) on a domain x Qi defining the contact between the 

qa 



sample and the l-th lead. The coupling strength tensor T^^ is modeled as 



= / drdv' {¥^{r')y *f (r)<?i^(r, r') + h.c. (7) 
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The integration domains for the leads are chosen to be 



= 

and for the system 
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with r G f^g and G defines the couphng of any two single-electron states by the 
'nonlocal overlap' of their wave functions in the contact region of the leads and the 
system, and their affinity in energy. A schematic view of the coupling is presented in 
Fig. [l] The parameters 5[ and 5[ define the spatial range of the coupling within the 
chosen domains f^g x Qi. We will be analyzing systems with complex subband structure 



H^t) 



Hjt) 



Figure 1. A schematic view of the coupling of the system to the leads. The green 
shaded areas correspond to the contact regions defined by the nonlocal overlap function 
^L,R in Hrit). 



and geometry. For that purpose we need quite many single electron states (SESs) that 
in turn lead to a requirement of an exponential number of many-electron states (MESs). 
In order to deal with this computational problem we define a window of relevant SESs 
with energies in the interval [fiR — A, //i, + A] that will consequently be used to build the 
MESs used in the calculation, see Fig. [2j Whether, the window of SESs is large enough 
in a specific calculation can only be decided by numerical experimention. Below this 
window all states are considered permanently occupied and above it all states remain 
empty (see Ref. [20]). 

As has been described in detail earlier ^20j we shall deploy an occupation number 
basis constructed from the SESs of the isolated finite quantum wire {^^^}- The MES // 
is then 

|/.) = |^^^^,...,^^:,...), (11) 
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Figure 2. A schematic view of the chemical potential bias of the leads and the relevant 
single-electron states of the system (Red). A/i is the bias window. 



where the number i^ is 1 if the n-th SES is occupied and if it is empty. We use here 
Greek letters for the labelling of the MESs and Latin letters for the labelling of the 
SESs. According to the selection of the relevant SESs around the bias window they will 
be of the form 

l/i) = I M^, 4o+i, 0, 0, ), (12) 

A^o states 

where A^o denotes the highest SES below the bias window, but A^max is the number of 
the highest SES inside it. 

The many-electron statistical operator of the whole system (leads and quantum 
wire) obeys the Liouville-von Neumann equation 

iWit) = [Hit), W{t)], Wit < to) = PlPrPs, (13) 

where the equilibrium density operator of the disconnected lead / with chemical potential 
ILii is 

^ Tr/{e-^(^^-^^^^)}' ^^^^ 

ps is the density operator of the isolated quantum wire {t < to) and serves as an initial 
condition for the reduced density operator (RDO) defined as the (partial) trace over the 
Fock space of the leads 

p{t) = TT^TT^W{t), p(to) = Ps. (15) 

We do not impose equilibrium condition on the leads after they have been coupled to 
the system, the finite quantum wire, at t > to- This is different from other approaches 
where one imposes W{t) = PLPRP{t)^ (see for example Ref. ||9]). In the second order in 
the coupling Hamiltonian the time evolution of the RDO is governed by the operator 
integro-differential equation [20j 

m = -l[Hs,pit)] - ^ E / dqx'it){[T',^At)] + h.c.), (16) 
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where we have introduced two operators to compactify the notation 



^qiit) = e-''''^ / dsx\s)Uqi{s)e'^'-'^'^''^e''''^, (17) 

J to 

U,i{s) = e-^« (T'tp(,)(l - /') - p(s)T't/) e-^^ (18) 
and a scattering operator T acting in the Fock space of the system 

a,f3 a 

T^p{q) describes the 'absorption' of electrons from the leads to the system and changes 
the many-electron state of the latter from f3 ^ a. The Fermi function of the SES 



labelled by g ^ i^yQ) l^^d / is noted by / . In order to derive Eq. (|16|) we have 



used the projection operator V = PR{to)pL{to)TrRTrL to project out an equation for 
the evolution of the system under the influence of the leads. Though we do not do 
it here, we could equally well project out an evolution equation for the leads. This 
equation would tell us how the electronic system in the leads evolves for t > to without 
assuming the central system to stay in equilibrium. Thus the Fermi distribution present 



in Eq. (18) is valid only for t < to and the eflFective distribution for later times can 
only be found from the RDO for the leads. This ansatz is not expected to change the 
steady-state currents, but it will influence the transient regime. 



Projecting (16) on the many-electron states constructed from the relevant single- 
electron states in the extended bias- window results in A/" = 2^™^^"^° coupled integro- 
differential equations for the matrix elements {a\p{t)\(3) of the reduced density operator 
describing the time evolution of the open system, i.e. the flnite quantum wire coupled 
to the leads. 

With the RDO it is now possible to calculate the statistical average of the charge 
operator Qs = e rf^rf^ in the coupled quantum wire 

{Qs{t)) = TT{W{t)Qs} = Trs{[TrLRVr(t)]gs} 

= Trs{p{t)Qs} = eJ2€ (Mm, (20) 

with the traces assumed over the Fock space. We are also interested in the average 
spatial distribution of the time-dependent charge 

{Qs{r,t)) =eJ2J2'^:ir)%{r)p,^{t){u\SM. (21) 

ab iJih> 

The net current flowing into the sample is 

HMt)) = {Jh{t)) - {j^{t)) = ^^^^ = eJ2Y.'aWim ■ (22) 

a jji 



Through the GME (16) it is possible to identify the contribution of each SES in the 



system to the current from the left lead or into the right lead |20j. More precisely, Eq. 



(16) gives us an expression for p{t). The trace of the flrst part, the commutator of p 
and i^s vanishes. The summation over leads in the second part allows us to identify the 
contribution to or from each lead. 
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The iterative scheme to solve the GME has been described in an earher pubhcation 



[20] . Due to the commutator structure of the GME (16) the conservation of probabihties 
can be verified, i.e. Trsp(t) = 0. At any time step in the numerical iteration we check 
the conservation of probability. 

3. Results and discussion 

An advantage of the GME formalism is the freedom to select the initial state of the 
system in any way compatible with the description via the RDO. We can thus choose 
the system to be initially empty, or occupied by any number of electrons infiuenced by 
the equilibrium state in either lead. More complex nonequilibrium initial states are also 
possible. 

In our calculation we assume GaAs parameters, m* = O.OGTmg, = 12.4, and 
set to = 0. The characteristic energy of the parabolic confinement for both the finite 
quantum wire and the leads is assumed HQq = 1.0 meV. The coupling between the leads 
and the quantum wire is described by the function 



X^'^(t) = <' ^-''"hJ (23) 

[ if t < 

with = 1.0 ps~^. We fix the temperature of the reserviors at T = 0.5 K. The 



parameters determining the coupling of the subsystems in the function g^^^ (10) are 
5ia^ = 1.0, 52<^^ = 2.0, and = 1.0 meV in order to let the coupling strength be 
determined by the behavior of the eigenstates of the quantum wire and the leads only 
very close to the contact region at x = ±Lx/2. For the numerical coupling constant g^ 
in (10) we select the value g^^aJi^ = g^^aJi^ with g\^ = 10.0 meV for both leads, or a 
different value stated for a particular calculation. (The rather unusual dimension of the 
numerical coupling constant comes from the fact that the x-part of the wave function in 
the leads is only (^-normalizable, see Appendix A). The numerical value of ^^qo P^^ 
se a limited significance, but below we shall explicitly show the effective coupling of the 
states in the leads and the system graphically. 

A time step At = 0.01519 ps was selected in order to allow for time integration 
to relatively large times without any visible loss of accuracy on the scales presented 
in the figures to follow. As will be shown in subsequent figures we include part of 3- 



4 energy subbands in the g-integration in the GME (16). For the wavevector cut-off 
qayj = 4.0 it has proven essential to use a large enough number of ga^-integration points 
to guarantee high accuracy for large times. Within each subband we use a 4-point 
Gauss integration repeated 180 times. This requirement weighs heaviest in increasing 
the needed CPU-time. 

Embedded in the finite quantum wire we have a quantum dot or an antidot 
represented by the Gaussian potential 

V{x) = Vo exp {-[(3,{x - x,)f - [(3y{y - yo)f]. (24) 

with = = 0.03 nm~^, and Vq and the spatial off-sets to be stated later. 
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3.1. Pure finite wire 

Before presenting results for the quantum wire with an embedded subsystem we use 
results for the time-dependent transport through the pure finite quantum wire to 
familiarize us with the system. 

The energy spectrum of the leads is shown in Fig. |3] together with the chemical 
potential in each lead, fiL = 0.9 meV, fiR = 0.7 meV, and the limits fiR — A and 
/LiL + A with A = 0.1 meV defining the window of relevant states around the applied 
bias e]4ias = I^l — f^R = 0.2 meV. The maximum energy for each subband shown in 
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Figure 3. The energy spectrum of the leads (solid red) vs. the scaled wavevector qa^^ 
and the chemical potential in each lead jj^l = 0.9 meV, jur = 0.7 meV, and the window 
of relevant states [fiR — A, /i/, + A] for A = 0.1 meV, HQq = 1.0 meV. 

the graph indicates the corresponding maximum wavevector in the ga^-integration of 



the GME (16). The energy spectrum of the 900 nm long quantum wire is shown in 
Fig. [ij together with the extended bias window [jj.R — A, fi^ + A] containing 5 SESs, 
the relevant SESs for including in the transport calculation and construction of the 
MESs. The actual bias window contains 2 SESs. The finite quantum wire is long 
enough compared to its effective width to show a clear indication of formation of energy 
subbands in the energy range shown in the figure. 

The time-dependent occupation of the relevant SESs is shown in Fig. [5] in 
comparison with the time-dependent coupling function x^'^- 

The numbering of the relevant SESs with the index a is in the order of increasing 
energy. In the initially empty system the many-electron state is labelled by /xq = 1, and 
the two-electron state with the electrons in the two lowest SES is labelled by /xq = 4 
[20j. We note that some of the higher SES seem to reach a steady state fractional 
filling fast, while the ones lower in energy are still increasing their occupation at times 
as large as 120 ps. Two effects contribute to this, the higher lying states are stronger 
coupled to states in the leads as we will show below, and they can conduct faster. The 
time-dependent GME formalism introduces energy dissipation from the finite quantum 
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Figure 4. The energy spectrum of the system (red dots) vs. the SES state number 
and the chemical potential in each lead /iL = 0.9 meV, jj^r = 0.7 meV, and the window 
of relevant states [jiiR — A, /ii, + A] for A = 0.1 meV. = 900 nm, HQq = 1.0 meV. 




Figure 5. The time dependent occupation of the relevant SESs for an initially empty 
system /io = 1 (a), and for a system occupied with initially 2 electrons {jno = 4) in 
equilibrium (b). The time coupling function is shown for reference. Lx = 900 nm, 
hfto = 1.0 meV, gl^Q = 10.0 meV, A^ = 1.0 meV, and T = 0.5 K. 

wire into the leads. At the finite temperature, T = 0.5 K, we see that in case of two 
electron initially in the system there is a finite but small probability for the electrons to 
get out of the system, even though they have energy below the actual bias window. As 
expected, the system looses the electron occupying the SES closer to the bias window 
with a higher probability. 

The time-dependent total current injected into the system from the left lead J^{t) 
and the total current leaving the system into the right lead are shown in Fig. 

[6| Here we see that in case of the initially empty system, /xq = 1, the current in the 
right lead, J^, is negative meaning that it is directed into the system for t < 40 ps. 
The system is supplied with electrons from both ends initially before it reaches a steady 
state with constant current through it. In case of the system initially occupied by two 
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Figure 6. The total current entering the system from the left lead J^{t)^ and the 
total current exiting the system into the right lead J^{t) vs. time for an initially 
empty system /io = 1 and a system with two electrons initially in equilibrium /io = 4. 

= 900 nm, hQo 1.0 meV, g^Q 10.0 meV, = 1.0 meV, and T 0.5 K. A 
positive sign of the current indicates a flow of electrons with charge e from left to right. 



electrons, /xq = 4, we do not have a net transfer of charge from the right lead initially, 
but we see fluctuations in for t < 20 ps before it turns positive when a net current is 
flowing through the system. Not surprisingly, the current has a maximum value when 
the occupation is changing the fastest for the system, see Fig. [5| 

The coupling strength tensor introduced earlier ^ gives the coupling between a 
state a in the relevant extended bias window and a state qUy in the leads (here the 
coupling is the same for left and right lead). In Fig. [7] we see the tensor for the 
3 lowest subbands of the leads labelled by Uy. The maxima correspond to resonant 
tunneling when Ea = e^^(g). The type of coupling selected ^ reproduces an eflFect 
well known in models for multimode transport in quasi-one dimensional quantum wires 
built on the Lippmann-Schwinger scattering approach [23j, i.e. the parity of the subband 
wavefunctions in case of a symmetric system forbids coupling of the nearest neighboring 
subbands. Therefore, the coupling of the relevant states here all lying in the flrst subband 
of the flnite quantum wire to the second subband of the leads is vanishing, as is seen 
in the center subflgure of Fig. [7| The coupling to the third subband is reduced by the 
exponential term in energy of Eq. g. 

The RDO can be used to calculate the average spatial charge distribution (21) of 
the MES in the flnite quantum wire at any instant of time. In Fig. [8] we show it soon 
after the coupling of the system to the leads, and again when the system is close to 
reaching a steady state. Just as we have seen in Fig. [6] of the total current in the right 
and left leads, initially the probability density increases in the empty flnite wire from 
both sides with the higher bias to the left supplying it faster there. The steady state 
attained in the end is clearly a mixed state with contribution from all of the available 
SES, and the coupling to the leads maintains a higher probability at the ends of the 
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Figure 7. The coupling between a state a in the relevant extended bias window and 
state quy in the left (and the right) lead for Uy = (top panel), Uy = 1 (lower left 
panel), and Uy = 2 (lower right). = 900 nm, H^Iq = 1.0 meV, ^qq = 10.0 meV, 
A^P = 1.0 meV. 




-15 -10 -5 5 10 15 -15 -10 -5 5 10 15 
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Figure 8. The average spatial charge distribution for the MES constructed from the 5 
relevant SESs in the extended bias window for t = 0.76 ps (left panel), and t = 121 ps 
(right panel). Note the huge difference in scale. The system is initially empty /j^q = 1. 
900 nm, = 1.0 meV, g^Q = 10.0 meV, = 1.0 meV, and T 0.5 K. 

finite wire. 

3.2. Finite wire with an embedded subsystem 

We now continue our exploration of the effects of the geometry of the system and 
leads on the description of the time-dependent transport by the GME formalism by 
introducing a Gaussian potential into the finite quantum wire. In Fig. [9] the energy 
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spectrum of the SES is shown for an off-centered Gauss weU with parameters Vq = —2.0 



meV, and = = 0.03 nm for the potential given by Eq. (24). For aU the cases 



of an off-set Gaussian potential we are only considering an off-set in the ^/-direction 
and we shall keep Xq = 0. In the figure we indicate the extended bias- window for the 
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Figure 9. The energy spectrum of the system (red x) vs. the SES state number i for 
the system with an embedded off-centered Gaussian well compared with the spectrum 
of a purely parabolic system (blue +), and the chemical potential in each lead jj^l = 1.4 
meV, jiiR = 1.2 meV, and the window of relevant states [juji — A, /i/, + A] for A = 0.1 
meV. = 900 nm, hfto 1-0 meV, Vq -2.0 meV, (3x = Py = 0.03 nm-\ yo 30 
nm, and xq = 0. 

relevant states by A = 0.1 meV, but in the calculations to follow we will often use 
a larger value for A to be stated in each case. We start selecting the location of the 
bias window just below and touching the second subband of the system. Below we 
shall discuss the character of the SESs, but first we investigate the current through 
the SESs of the system. Experimentally only the total current entering or leaving the 
system is measureable, but the current through each SES can give us an insight into the 
transport processes in the system. The current from the left lead into each relevant SES 



is displayed in Fig. [10] for two different off-sets of the Gaussian well, the upper panel for 
the smaller off-set, z/q = 30 nm, and the lower panel for the larger one, i/q = 60 nm. We 
see that the current into the lowest state in energy, a = 1, which is below the chemical 
potential in both leads is highest the first 60 ps while it is reaching its steady state 
value. After this period the current reaches a steady state value that can be verified by 
checking that the current leaving the system into the right lead (not shown here) has 
the same value. The same can be said about the current through the states a = 3, 4, but 
the current through state a = 2 situated just below /llr does not seem to reach a steady 
state value in the time interval shown. Indeed, the current through a = 2, ^2, seems 
to be vanishing here. We shall investigate this further just below, but first we show in 



Fig. 11_ how indeed, the total current into the system from the left and the total current 
leaving the system on the right reach the same values soon after t = 60 ps. Here is also 
clear that during the first 20 ps the total current in the right lead is directed into the 




Figure 10. The time-dependent current from the left lead into the SES a for the 
system with an embedded off-centered Gaussian well at yo = 30 nm (left panel), and 
yo = 60 nm (right panel). A = 0.15 meV. = 900 nm, hQo = 1-0 meV, Vq = -2.0 
meV, Px = Py = 0.03 nm~^, ^qq = 7.5 meV, and xq = 0. 
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Figure 11. The total left and right current vs. time for the system with an embedded 
off-centered Gaussian well. A = 0.15 meV. = 900 nm, hQo = 1.0 meV, Vb = —2.0 
meV, Px — Py = 0-03 nm~^, ^qq = 7.5 meV, and = 0. A positive sign of the current 
indicates a flow of electrons with charge e from left to right. 



system supplying it with electrons. Though we do not show it here, the same is even 
true for all the partial right currents in this time interval, though strongly decreasing 
for SESs with higher energy. 

In order to investigate further the behavior of the current j2 through the second 
relevant SES in the extended bias window (see Fig. 10) we repeat the calculation for a 
stronger coupling ^^qo = 10-0 i^i^V, and a slightly wider extension of the bias window by 
using A = 0.2 meV. This change of the window only brings in two extra SES at the top 
of it. The results are seen in Fig. 12, where we now display both the left and the right 
partial currents. Like expected the occupation of the system here takes a shorter time 
and the currents are higher than for the case of lower coupling. There is of course a slight 
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Figure 12. The current from the left lead into the SES a (left panel), and the current 
from the SES a into the right lead (right panel) vs. time for the system with an 
embedded off-centered Gaussian well, yo = 60 nm, A = 0.2 meV. = 900 nm, 
hQo 1.0 meV, Vq = -2.0 meV, (3:, = Py = 0.03 nm-\ g^^ = 10.0 meV, and xq 0. 

rearrangement of the individual partial currents, but an eye catching change is that now 
the steady state values for and jf are reversed in comparison with the other partial 
currents. The partial current through a = 2 flows from right to left, but the total steady 
state current is still in the expected direction, from left to right. To get an idea why this 
is happening we present in Fig. [13] the probabilities of the six relevant SESs. There we 




Figure 13. The probability density of the single-electron eigenstates of the system 
labelled by a in numerical order with a = 1 at the top left along rows to a = 6 at the 
bottom right. In the system is an embedded off-centered Gaussian well. = 900 nm, 
hQo = 1.0 meV, Vq = —2.0 meV, yo = 60 nm, /3x = Py = 0.03 nm~^, and xq = 0. 

see that the SES a = 2 is a quasi-bound p-state in the off-centered Gauss well with thus 
a high probability in the well and just below it, but a much reduced probability density 
towards the contact ends of the system. In the well region it has a character of a state of 
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the second subband, but close to the contacts it has a character of a state from the first 
subband. In our earher calculations with the Lippmann-Schwinger model of wells in a 
totally open quantum wire such quasi-bounds states would always cause a sharp dip in 
the conductance due to a strong backscattering [23l EU [25]. Here, the oflF-centering of 
the Gaussian well is important in order to enhance the coupling between neighboring 
subbands. This explains the fact that the total current thought the system is slightly 



lower for the case of the more off-set well, (see Fig. 12). In calculations with the present 



GME model we have not found corresponding behavior for an embedded hill in the 
system in the case of no external magnetic field. The conclusion is thus close at hand 
that the reversal of the partial current here is a manifestation of the total refiectance of 
a quasi-bound state in the GME formalism. One might suspect the fact that the SES 
a = 2 is just below /xr playes a role here, but we have excluded this explanation by 
shifting the chemical potential window slightly down to place the state into its center. 
After this shift the same state still exhibits a reversed steady state current, so the 
character of its wave function plays the main role here and not its location in the energy 
spectrum with respect to /xl and /xr. 

Another surprise is in store when we look at the time-dependent average spatial 
charge distribution for the system in the beginning and after the system has reached 
a steady state displayed in Fig. [T4| At a very early time, t = 0.76 ps, we see electron 
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Figure 14. The average spatial charge distribution for the MES constructed from the 
6 relevant SESs in the extended bias window for t = 0.76 ps (left panel), and t = 68 
ps (right panel) for the system with an embedded off-centered Gaussian well. Note 
the huge difference in scale. The system is initially empty /io = 1- = 900 nm, 
hVtQ 1.0 meV, 10-0 meV, 1.0 meV, and A = 0.2 meV Vb -2.0 meV, 

7/0 = 60 nm, = 0, and Px = Py = 0-03 nm~^. 



probability seeping in from both contact regions, though more from the higher bias 
region at the left, and we see already probability entering the well. In the steady-state 
regime the system has already entered at t = 68 ps we have a strong p-state around the 
well, but the electrons have a higher probability to be found just outside the off-centered 
well. This may not be very surprising in light of the fact that the system carries a good 
amount of current through it, and we are looking at a MES here that both carries the 
information of electrons being quasi-bound in the well, and being scattered by the well. 
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This latter fact together with the bias difference between the ends causes the 'p-state' to 
be shghtly rotated from the y-axis. A further support for this picture can be sought in 
the results for an off-centered Gaussian hill shown in Fig. 15 for the same two moments 



of time explored for the case of the well in Fig. 14 



0.05 





Figure 15. The average spatial charge distribution for the MES constructed from the 
6 relevant SESs in the extended bias window for t — 0.76 ps (left panel), and t = 68 
ps (right panel) for the system with an embedded off-centered Gaussian hill. Note the 
huge difference in scale. The system is initially empty /io = 1- = 900 nm, Wt^ = 1.0 
meV, 7.5.0 meV, 1.0 meV, and A 0.2 meV Vb 2.0 meV, 30 nm, 

xo = 0, and Px = Py = 0.03 nm~^. 



Initially at t = 0.76 ps we see the charge seeping into the system from both contact 
regions as for the case of the well, but now no extra probability is seen close to the hill, 
it is a repulsive potential. When the system has reached a steady state at t = 68 ps we 
find out that the electrons have the highest probability to be found in the system close 
to the hill. The Gaussian hill is off-centered, thus mixing up the motion of the electrons 
along and perpendicular to the finite quantum wire, the system. Again, the average 
charge density is slightly tilted due to the external bias. Classically speaking we would 
say that electrons are scattered between the hill and the parabolic boundary of the wire 
opposite to the hill on the side of the wire with negative ^-coordinate. The electrons 
have to go through this symmetry breaking constriction and spend more time there 
than elsewhere in the system. It is interesting to view the probability distribution of the 
relevant SESs in Fig. [16] to see that this fact is even found in the stationary eigenstates 
of the system. In addition, we see that states a = 4, 5 with a character mainly refiecting 
the second subband and thus with higher energies suggest a classical analogue of zig-zag 
motion, and even the state a = 6 with a clear character of a SES in the first subband 
displays this zig-zag motion. A peek back at Fig. [13] reminds us that this behavior was 
already present for the system with an embedded well, though not quite as prominent. 

We have now seen that geometrical properties of the system affecting the 
bandstructure in the neighborhood of the bottom of the second subband, where, for 
example, a quasi-bound state can be found for the case of an embedded well. What 
about the bottom of the first subband? There the states are less coupled to the states 
in the leads and carry a smaller amount of a current. In Fig. [17] we show the energy 
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Figure 16. The probability density of the single-electron eigenstates of the system 
labelled by a in numerical order with a = 1 at the top left along rows to a = 6 at the 
bottom right. In the system is an embedded off-centered Gaussian hill. = 900 nm, 
hfto = 1.0 meV, Vb = 2.0 meV, yo = 30 nm, xq = 0, and Px = Py = 0.03 nm~^. 



spectrum of a system with a centered shallow well and the chemical potentials selected 
such that only the bound state of the well is below the actual bias window. 
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Figure 17. The energy spectrum of the system (red x) vs. the SES state number i 
for the system with an embedded centered Gaussian well compared with the spectrum 
of a pure system (blue +), and the chemical potential in each lead /iL = 0.55 meV, 
fiR = 0.45 meV, and the window of relevant states [/j^r — A, /j^l -\- A] for A = 0.136 
meV. 900 nm, = 1.0 meV, Vb -0.4 meV, and (3^ = f3y = 0.03 nm-^ 

We do not show the time-dependent occupation of the SESs here, but it would 
reveal the fact that even after the 121 ps the system is far from reaching a steady state, 
in fact the occupation of the levels with a > 1 is still growing linearly. This can also 
be verified by observing the partial left and right current for the relevant SESs in Fig. 



18, The current through the bound state ji remains negligible all the time and through 
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Figure 18. The current from the left lead into the SES a (left panel), and the 
current from the SES a into the right lead (right panel) vs. time for the system with 
an embedded centered Gaussian well, yo = nm, A = 0.136 meV. = 900 nm, 
hQo 1.0 meV, Vq = -0.4 meV, (3:, = Py = 0.03 nm-\ and 7.5 meV. 



the very small value of the probability for the bound SES a = 1 in the contact region a 
change in its occupation is not expected until on the nanosecond scale. Moreover, the 
system is here in the phase that some of the partial current in the right lead are still 
directed toward the system, supplying it with electrons. 



The time-dependent average spatial charge distribution is shown in Fig. 19 for 





Figure 19. The average spatial charge distribution at t = 121 ps for the MES 
constructed from the 5 relevant SESs in the extended bias window for initially empty 
system /io = 1 (left panel), and initially occupied by one electron fio = 2 (right panel) 
for the system with an embedded centered Gaussian well. = 900 nm, HQq = 1.0 
meV, 7.5 meV, 1.0 meV, and A 0.136 meV Vq -0.4 meV, = 

nm, and Px = Py = 0-03 nm~^. 

t = 121 ps indicating that the empty system is almost still empty at this time, or in 
the case of one electron initially in the system it is still there at this time without any 
significant change. The coupling between the bound state and the other states is simply 
much to small as it is governed by the behavior of the states in the contact region at the 
ends of the finite quantum wire. In contrast, the quasi-bound state corresponding to 
a = 2 in Fig. M (the second state from below in the window of relevant states marked by 
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the green dashed hnes) had a significantly higher couphng to other relevant transport 
states as can be seen in Fig. [13] from the fact that it has a higher probability in the 



contact region. It, was indeed a quasi-bound state, but the SES a = 1 shown in Fig. 19 



is a real bound state on the time scale shown (121 ps) with no discernible coupling to 
the relevant SESs above it in energy. The current through this state oscillates around 
zero with amplitude many orders of magnitude lower than the currents through the 



neighboring states seen in Fig. 18 



4. Summary 

In this publication we have shown that the GME formalism can be used to describe time- 
dependent transport through a semiconductor system on the nanoscale with complex 
geometry connected to broad leads with several active modes. We have focused our 
attention on internal processes in connection with bound and quasi-bound states in the 
system, rather than attempting to describe specific systems of experimental interest. 
In this way we have been able to compare the results of the GME formalism to earlier 
calculations with the Lippmann-Schwinger formalism j23l [2H [25] . The model presented 
here is very fiexible, in the sense that different parabolic confinement can be chosen for 
each part, the left and the right leads, and the central system. The coupling to each lead 
can be configured individually in time and space. In addition, the embedded system or 
even the shape of the central system can be changed without a problem. All this gives 
interesting possibilities that have to be explored in the future work in closer contact 
with experimental work. 

Geometrical eflFects like selection rules for processes between the subbands of the 
system and the leads have been successfully incorporated in the model by selecting a 
coupling of the type represented by Eq. ([t]). We have found the resulting dynamic 
current through the system to depend quantitatively on the "size" of the contact area 
by varying S[ and but the qualitative behavior of the current is not very sensitive to 
variation of these parameters within the same magntitude. 

It is only fair to suggest possible experimental systems or setups at this stage where 
our GME model could be tested. The first direct implimentation could be a long broad 
quantum wire with the coupling region, indicated by green or shaded areas in Fig. [T| 
realized by gate strips situated above the lead-sample contact regions, controlled by 
external potentials. Here, like in our model the subband structure in the leads and 
the system would have a large infiuence on the transport through the system. Chen et 
al. [i26j have measured the static magnetotransport properties of a quasi- ID quantum 
wire where one might investigate whether the coupling to the broad ,, leads" is better 
described by our GME model appropriate for weak coupling or the Lippmann-Schwinger 
approach appropriate for a more ,,open" access or coupling to the leads. Definitely, 
a more developed GME model could be appropriate to explore the geometrical and 
dynamic properties of the single quantum dot studied by Astafiev et al [27j. 

The numerical accuracy of our results presented here has been tested by including 
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more SESs in the calculation of the states and the energy spectrum of the system, and 
by including more subbands in the leads and denser g-points in the integration of the 
GME. We have carefully selected the time-mesh fine enough for the GME and where 
possible we have experimented with different sizes for the window of relevant states 
determined by A. 

We are aware that the time-dependent properties of the coupling have to be selected 
in accordance with the relatively narrow energy range of the MESs constructed for the 
central system. This set of MESs may only be appropriate for time-dependence that is 
not too fast or strong with respect to the window of relevant SESs. 

Moreover, it is clear that electron-electron interaction effects neglected here may 
be of paramount interest in experimentally relevant systems. Our only excuse is that 
we have here taken the first steps to use the GME formalism for a system with rich 
geometry without resorting to the Markov approximation. Inclusion of the Coulomb 
interaction is in no way trivial. We are dealing with a system with a variable number of 
electrons where different charging regimes may be of importance depending on the type 
of coupling between the leads and the system. The huge number of MESs needed for 
the description of the Coulomb interaction in order to retain correlation effects (that is, 
going beyond a mean field approximation) is a real problem and opens the important 
issue of how to isolate a relevant section of the Fock-space. This is a nontrivial future 
task that we will not shy from. 

It is well known that the GME used here does not guarantee the reduced density 
matrix to be positive definite unless the coupling between then leads and the central 
system is sufl&ciently weak. We have checked this for each iteration step in our 
calculation and show only results here for a positive definite density matrix. For states 
high in the subbands negative probabilities may occur earlier than for states close to 
the bottom of the first subband, as this phenomenon depends on the effective coupling 
of the relevant states to states in the leads. 
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Appendix A. The single-electron states of the system and the leads 

The SESs of the finite quantum wire and their energy spectrum are found by 
diagonahzing hs ((T|) in the basis {(Pn^<j)ny} with 



^cosf™V n^ = 1,3,5,... 

^nM=< , (A.1) 

satisfying the hard wah boundary condition at x = ±Lx/2 and 

2a f I, / y 



The matrix elements of the embedded potential (24) have been calculated analytically. 

For the semi-infinite leads with the same confinement parameters we can use the 
same basis functions for the ^-direction, but for the x-direction we use 

^q{x) = sin [q{x ± (A.3) 

V ZTT 



with in the left lead and ' — ' in the right lead, {^fq} is a complete orthogonal 
and 5-normalizable basis for the continuous wavenumber q > 0. These eigenfunctions 
represent the fact that in equilibrium before the system and the leads are coupled 
together the states along the leads are standing waves with an equal amount of left and 
right propagating waves. 
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